A systematic comparison of jet quenching in different fluid-dynamical models 
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Comparing four different (ideal and viscous) hydrodynamic models for the evolution of the medium 
created in 200 A GeV Au-Au collisions, combined with two different models for the path length de- 
pendence of parton energy loss, we study the effects of jet quenching on the emission-angle depen- 
dence of the nuclear suppression factor Raa{4>) and the away-side per trigger yield Iaa{4>)- Each 
hydrodynamic model was tuned to provide a reasonable description of the single-particle transverse 
momentum spectra for all collision centralities, and the energy loss models were adjusted to yield the 
same pion nuclear suppression factor in central Au-Au collisions. We find that the experimentally 
measured in- plane vs. out-of-plane spread in _Raa(0) is better reproduced by models that shift the 
weight of the parton energy loss to later times along its path. Among the models studied here, this 
is best achieved by energy loss models that suppress energy loss at early times, combined with hy- 
drodynamic models that delay the dilution of the medium density due to hydrodynamic expansion 
by viscous heating. We were unable to identify a clear tomographic benefit of a measurement of 
Iaa{(I>) over that of i?A,4(0). 

PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 

The expression 'jet tomography' is often used to de- 
scribe the analysis of hard pQCD processes taking place 
inside the soft medium created in an ultrarelativistic 
heavy-ion collision, with the aim to study properties of 
the medium. In particular, the focus is often on the nu- 
clear suppression of hard hadrons in A-A collisions com- 
pared with the scaled expectation from p-p collisions, due 
to loss of energy from the hard parton by interactions 
with the soft medium (see e.g. [IHa]), expressed through 
the nuclear suppression factor Raa- 

In comparing theoretical calculations with experimen- 
tal data on Raa, there are two main unknown properties 
of the medium: The nature of the parton-medium inter- 
action, being closely connected with microscopic proper- 
ties of the medium (such as the relevant degrees of free- 
dom) , and the evolution of the medium density distribu- 
tion, being connected with macroscopic properties such 
as the thermodynamical parameters in a fluid descrip- 
tion of the medium. While some attempts at systematic 
comparison of different models for the parton-medium 
interactions using the same fluid-dynamical model for 
the medium have been made in order to assess the effect 
of assumptions in the parton-medium interaction model 
[1, 0: there is very little systematies available for the 
effect of different hydrodynamical models on jet quench- 
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ing observables other than the overall suppression ratio 

RAAipx) UM- 

What may have slowed the insight that there is a need 
to systematically understand the role of the medium den- 
sity evolution is the fact that early comparisons with data 
were usually done on the basis of single-hadron suppres- 
sion Raa for central collisions only, and it took some time 
before it was realized that this quantity is quite insensi- 
tive to model assumptions [l3l, especially when (as usu- 
ally done) one model parameter governing the strength 
of the parton-medium interaction is fit to the data. The 
need for more differential observables, such as RAAift^) 
as a function of the angle of the observed hadron with 
the reaction plane for different centralities [ll| or the 
strength suppression Iaa observed in hard back-to-back 
correlations [8|, to overcome this insensitivity was only 
realized later. 

Such observables are primarily sensitive to the effec- 
tive path length dependence of the energy loss. As 
one goes from central to peripheral collisions, both the 
mean density of the medium and the average path length 
needed for a hard parton to traverse the medium de- 
crease. Within a given centrality class, Raa{4') is domi- 
nated by the change in path length, modulated by a weak 
directional dependence of the average density probed by 
the parton. How precisely the path length and density 
change with centrality depends, however, on details of 
the hydrodynamical evolution. 

The aim of this paper is to investigate in some de- 
tail the connection between high-py observables and the 
bulk medium evolution. In particular, we try to iden- 
tify those properties of a hydrodynamical model which 
have the strongest influence on high-pT observables. We 
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do so by presenting a systematic study of the direc- 
tional dependence of the nuclear suppression factor Raa 
and the away-side yield in triggered back-to-back correla- 
tions, Iaa, for several parton-medium interaction models 
with different path length dependence and a number of 
different hydrodynamical models for the medium. The 
hope is to derive constraints for a combination of both 
medium evolution and parton-medium interaction mod- 
els that can be used to eventually arrive at a detailed 
understanding of the dynamics of ultrarelativistic heavy- 
ion collisions. 



II. HYDRODYNAMICAL MODELS 

We describe the medium probed by the hard parton 
as a thermalized fluid. Its temperature and energy and 
particle densities evolve in space and time, due to hy- 
drodynamic expansion driven by pressure gradients. In 
this work we use both ideal and viscous hydrodynamics 
to generate these density profiles. 



A. Ideal hydrodynamics 



In the ideal case we solve the hydrodynamic equations 



(1) 



where T*^" = (e + P)u^u'^ — g'^'^P is the stress-energy 
tensor, jg — ubu^ is the baryon number current, ub is 
the net baryon number density, e the energy density, and 
P the pressure in the local rest frame which moves with 
fluid four-velocity in the global frame. The Equation 
of State (EoS) P = P(e, n^), relating the pressure to the 
local energy and net baryon number density, closes the 
set of dynamical equations. 

For testing parton energy loss, we have at our disposal 
space-time profiles of e, P and temperature T from two 
different ideal hydrodynamical models. The first of these 
(l^ solves Eqs. ^ in 3-t-l dimensions, propagating both 
T^^ and jg- The second model [l^l simplifies the prob- 
lem to 24-1 dimensions by assuming longitudinal boost 
invariance (i.e. none of the physical quantities depend 
on space-time rapidity = ^ h\[{t+ z) / {t— z)]) and setting 
the net baryon density everywhere to zero (such that only 
the energy-momentum tensor T^^ needs to be evolved, 
using a simplified form P(e) for the EoS). These approx- 
imations can be made since we are interested in energy 
loss only at mid-rapidity where, at RHIC energies, the 
net baryon density is very small. Both calculations use 
light-cone coordinates (r, x, y,ry), where r = \fW—z^ is 
the longitudinal proper time and t] is the space-time ra- 
pidity. 

Both models use smooth energy density distributions 
as initial conditions, based on the densities of binary col- 
lisions and wounded nucleons For details we re- 
fer to the original papers describing the models [H, [l^ . 



The hydrodynamic evolution starts at initial time tq = 
0.6 (0.17) fm/c in the (3-hl)-d ((2-hl)-d) model. The 
(3+l)-d model uses a bag model EoS with a first or- 
der phase transition at = 160 MeV [l^l whereas the 
(2-hl)-d model [11] uses the EoS from Ref. [Hi. Both 
Equations of State assume chemical equilibrium among 
the hadrons in the dilute resonance gas phase below T^. 
Thermal hadron spectra are calculated using the conven- 
tional Cooper-Frye method (l6j . where particle emission 
is calculated from a constant-temperature surface. The 
freeze-out temperature is Tdoc = 130 (160) MeV for the 
(3-|-l)-d ((2-|-l)-d) model. Strong and electromagnetic 
two- and three-particle decays of unstable hadrons are 
taken into account before comparing with experimental 
data. 



B. Viscous hydrodynamics 

We also study parton energy loss in a medium whose 
space-time evolution is computed from viscous hydrody- 
namics, by solving the second-order Israel- Stewart equa- 
tions in 2-1-1 dimensions as described in Ref. [l7[, assum- 
ing longitudinal boost invariance and zero net baryon 
density. Here the energy-momentum tensor of the fluid 
is decomposed as 



(2) 



which differs from the ideal fluid decomposition in 
Sec. Ill Al by the appearance of the traceless and symmet- 
ric shear viscous pressure tensor tt''^ satisfying u^tt^'' = 
0. Effects from bulk viscosity are neglected as small com- 
pared to ■n^'^ [ist . The energy- momentum conservation 
equations d^T^^ ~ arc supplemented by the Israel- 
Stewart [l^, [lOl evolution equations for the viscous pres- 



sure components tt'''^, see [17[ for details. 

The viscous hydrodynamic energy density profiles 
studied here were obtained with the Equation of State 
s95p-PCE described in [2l[ and Appendix C of Ref. [12]. 
It matches the latest lattice QCD data of the EoS at 
high temperatures with a chemically frozen hadron reso- 
nance gas EoS at low temperatur es l2ll . l2^ that uses non- 
equilibrium chemical potentials [23| to ensure preserva- 
tion of the stable hadron ratios at their chemical freeze- 
out values as the system cools below the chemical de- 
coupling temperature Tchom = 165 McV that has been 
experimentally established [1^ . 

In the viscous simulations we start the hydrodynamic 
evolution at tq = 0.4 fm/c and decouple the hadron mo- 
mentum spectra at Tdoc = 130 MeV. To compute the 
hadron spectra from the hydrodynamical output along 
the freeze-out surface we again use the Cooper-Frye pre- 
scription, but with a modified expression for the distribu- 
tion function, f{x,p) = fcq{x,p) + Sf{x,p), where we add 
to the local equilibrium distribution a small viscous cor- 
rection Sf that depends on the viscous pressure compo- 
nents 'k'^'^{x) at freeze-out and grows quadratically with 
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Pt (see [13, HH for details) . The specific shear viscosity 
is fixed at rj/s = 0.2, independent of temperature. 

In addition to the Glauber model initial conditions 
used in the ideal fluid dynamical models, we also study 
a set of viscous hydrodynamic evolution models based 
on Color Glass (CGC-fKLN) initial conditions [M 
which, for noncentral collisions, feature somewhat larger 
initial eccentricities and surface density gradients than 
the Glauber model profiles (see Fig. 1 in [2l[ for a com- 
parison of these profiles) and thus generate more ra- 
dial and elliptic flow. We will label viscous hydrody- 
namic simulations initiated with Glauber model profiles 
as "vGlb", and those initiated with CGC-fKLN profiles 
as "vCGC". We will discuss the consequences of these 
differences on the directional dependence of energy loss 
suffered by a parton propagating through these fireballs. 

Neither the ideal nor the viscous fluid simulations stud- 
ied here account for event-by-event fluctuations of the ini- 
tial shape and orientation of the collision flreball. Such 
calculations were recently reported both for ideal p?! - 
[29j and viscous hydrodynamics [soj . Source eccentrici- 
ties and anisotropic flow are most strongly affected by 
these fluctuations at very small and very large impact 
parameters where the flreball is either almost round or 
very small. Neither of these two situations is of interest 
in our present study. 



C. Spectra and elliptic flow from the 
hydrodynamic models 

Whereas the directional dependence of the soft hadron 
spectra (i.e. their elliptic flow) reflects the momentum 
anisotropy of the hadron emitting source at freeze-out, 
the emission angle dependence of parton energy loss 
probes more directly the geometrical aspects of the fire- 
ball, i.e. its spatial deformation. Still, it may matter 
whether the hard parton moves with or against the col- 
lective flow as it propagates through the fireball, so the 
hydrodynamical models need to be tuned to give a rea- 
sonably accurate representation of the momentum-space 
structure of the fireball, as reflected in the flnal hadron 
spectra, before we test the influence of differences in their 
geometrical features on parton energy loss. 

In this subsection we demonstrate that all models pro- 
vide a reasonable description of the transverse momen- 
tum spectra of pions and protons over the entire range 
of collision centralities jgO] . Testing the spectra for both 
a very light and a heavy hadron species, which react dif- 
ferently to radial collective flow [3l|, ensures that the 
appropriate amount of radial flow is generated in the evo- 
lution. The elliptic flow coefficient tests whether we also 
have the correct amount of flow anisotropy. Expressed 
in terms of velocity differences, the flow anisotropy is a 
small effect superimposed on a much larger radial flow 
velocity; while this anisotropy is crucial in determining 
the transport coefficients (in particular the shear viscos- 
ity) of the fireball fluid, it is not expected to lead to major 



modifications of the directional dependence of parton en- 
ergy loss. It affects the latter mostly by influencing the 
evolution of the spatial deformation of the flreball. As we 
will see, the time dependences of the radius and shape 
of the fireball play minor roles in the parton energy loss; 
hence, it is not a serious problem for our analysis that 
the elliptic flow V2(j>t) is not well reproduced by some of 
the hydrodynamical models we have studied. 

Figure [T] shows the pion and proton spectra for 
200 AGeV Au-Au collisions at various centralities, as ob- 
tained from the three different (2-|-l)-d hydrodynamical 
models. Except for the peripheral bins, the viscous sim- 
ulations give slightly steeper pion pT-distributions than 
the ideal fluid ones. This is a consequence of the differ- 
ent chemical composition in the hadronic phase as de- 
scribed in the preceding subsection [2^ [s^l . The proton 
spectra arc markedly better described by viscous than 
by ideal hydrodynamics. Shear viscosity adds to the ef- 
fective transverse pressure, generating more radial flow 
(TtI and thus making the spectrum harder, as desired by 
the experimental data (6]| . The somewhat steeper sur- 
face density gradients of the initial CGC energy density 
profile when compared to the Glauber model adds an 
additional small contribution to this effect. 

The pT-dependent pion and proton elliptic flow from 
minimum bias Au-Au collisions at y/s = 200 A GcV is 
shown in Fig. [2] Here the viscous calculations are seen 
to badly underpredict the measured values, as a conse- 
quence of shear viscous suppression of fiow anisotropics 
[17. . 35, 36]. The disagreement is not quite as bad for 
the CGC initial conditions which have somewhat larger 

initial spatial eccentricity Sx = |^2^i^2| [13, llll where 
((. . . )) denotes an average over the energy density in the 
transverse plane. But even in this case, the chosen value 
T]/s = 0.2 yields too much suppression of V2, in agree- 
ment with the findings in |37l |. We found in [2l| that 
quite generally the measured fiatness of the proton pT- 
spectra and the large elliptic fiow V2{pt) are in tension 
with each other, the former preferring larger rj/s whereas 
the latter likes smaller viscosities. We here chose to opti- 
mize the slope of the proton spectra, i.e. the magnitude 
of radial fiow. 



D. Geometry of the hydrodynamic fireballs 

The different initial and final conditions, equations of 
state and viscosities assumed in the different hydrody- 
namical fireball models lead to some differences in the 
space-time evolution of the size and shape of the fire- 
ball which, weighted with the corresponding density dis- 
tributions, affect the parton energy loss and its direc- 
tional dependence. Figure |3] shows the freczc-out sur- 
faces for Au-Au collisions in the 20-30% centrality class 
(b = 7.49 fm) in the r-r and x-y planes. The different 
shapes of the freeze-out contours along the x and y direc- 
tions in Fig.[3Ja) reflect the elliptical source deformation 
in non-central collisions. Figure EJ^b) shows freczc-out 
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FIG. 1: (Color online) Transverse momentum spectra of positively charged pions (a) and protons (b) from ideal and viscous 
(2+l)-d models for 200 A GeV Au-Au collisions at different centralities. Data from the PHENIX Collaboration [s^ are shown 
without error bars since errors are smaller than the symbol size. 



0.2 



0.1 



0.0 



2+ld ideal 


1 1 1 1 1 1 1 1 1 1 1 1 1 


2+ld vCGC 


(a) 


2+ld vGlb 


□ PHENIX 




TT 




□ 

D / 


. — _ . - ' ' 




Au+Au 




200 A GeV 


□/-<-' 

- V''' 

1 1 1 1 1 1 1 1 1 1 1 


Minimum bias _ 
1 1 1 1 1 1 1 1 1 1 1 1 1 



0.0 0.5 1.0 1.5 2.0 

Pt [GeV] 



2.5 



3.0 



0.2 



0.1 



0.0 



2+ld ideal 

2+ld vCGC 

2+ldvGlb 

□ PHENIX 


1 1 1 1 1 1 1 1 1 1 1 1 1 

(b) 


: p 








/ y"^ ' ' ' ' 

/ ° Au+Au 
/ n°o^-'' 200 A GeV 
/ ^ Minimum bias _ 


-utPr'T'l 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 



0.0 0.5 1.0 1.5 2.0 

Pt [GeV] 



2.5 



3.0 



FIG. 2: (Color online) Elliptic flow of (a) charged pions and (b) protons in minimun bias 200 A GeV Au-Au collisions. Data 
from the PHENIX Collaboration [i^] are shown without error bars since errors are smaller than the symbol size. 



contours in the x-y plane at times r = 2fm/c (solid) and 
r = 4fm/c (dashed), illustrating how the out-of-plane 
elongation of the source decreases with time [6^ • 

The most striking feature of these contour plots is 
the similarity of the freeze-out contours for the three 



(2+l)-dimensional models. Even though the starting 
time To for the viscous simulations (0.4fm/c) is more 
than twice as large as that used in the (2+l)-d ideal runs 
(0.17 fm/c), and the viscous fluid is allowed to cool down 
to Tdec = 130 MeV (compared to Tjec = 160 MeV for 
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FIG. 3: (Color online) Freeze-out surfaces of the different hydrodynamic models in the (a) r-r and (b) x-y planes for 200 A GeV 
Au-Au collisions in the 20-30% centrality class. The two halves of panel (a) show cuts through the freeze-out surface along 
the in-plane [x) and out-of-plane (y) directions, respectively. In panel (b), solid (dashed) lines represent cuts at r = 2fm/c 
(4fm/c). The contours for the three (2-|-l)-d hydrodynamic models are almost indistinguishable. 



the (2-|-l)-d ideal runs), all three models complete their 
freeze-out almost at the same time 7.2 — 7.5fm/c for 
h = 7.49 fm). This is due to the additional radial flow 
produced by shear viscous pressure and the use in the 
viscous flow simulations of EoS s95p-PCE vifhich is stiffer 
around the phase transition than the EoS used in the 
ideal fluid simulations [2l|. The combination of these 
two effects allows the viscous fireball to cool faster dur- 
ing the late evolution stages and thus freeze out sooner 
than an ideal fluid with the same initial conditions [TtI . 

During the early stages, on the other hand, viscous 
heating delays the cooling process, so at early times (up 
to about 4fm/c) the fireball center remains hotter and 
denser in the viscous case than for the ideal fluid [63| . 
The somewhat steeper initial density gradients of the 
CGC-fKLN profile (see Fig. 1 in (U) generate slightly 
larger radial flow in the vCGC model, causing it to fully 
decouple a fraction of a fm/c earlier than the vGlb model. 

Compared to the (2-|-l)-d simulations, the space-time 
volume covered by the ideal (3-|-l)-d fluid is much larger 
(wiggly green lines in Fig. This not a consequence 
of dramatically different transverse expansion in (3+l)-d 
and boost-invariant (2+l)-d evolution (near midrapidity 
the longitudinal density profiles and expansion velocities 
in the (2-t-l)-d and (3+l)-d simulations are very similar), 
nor is it primarily due to starting the (3-|-l)-d simulation 
later (at tq = 0.6 fm/c). The main reason is that the 
(3-|-l)-d simulations use a bag model EoS with a first 
order phase transition that produces a relatively long- 
lived mixed phase where the speed of sound vanishes and 



the fluid stops accelerating. 



III. PARTON-MEDIUM INTERACTION 
MODELS 

In this work, we use two different models for the 
parton-medium interaction. The first is the Armesto- 
Salgado- Wiedemann (ASW) model of medium-induced 
radiative energy loss in perturbativc Quantum Chro- 
modynamics (pQCD), in the formulation of energy- loss 
probability distributions, so-called 'quenching weights' 
|41| . As characteristic for perturbativc models of 
medium-induced radiation, it leads to a quadratic de- 
pendence of mean energy loss with the in-medium path 
length L in a constant medium due to the LPM suppres- 
sion of subsequent radiation processes. We have picked 
this model among other formulations of radiative energy 
loss, since it shows the strongest path length dependence 
for RAA{<f>) of all radiative energy loss models tested in 
the same hydrodynamical background Q • 

The second model is based on strong-coupling ideas 
for the medium. It is a hybrid model, in which the hard 
scales in the process are treated perturbatively, as in the 
standard pQCD radiative energy loss calculations, while 
the interaction with the plasma which involves strong- 
coupling dynamics is modeled by AdS/CFT calculations 
for the N = A super- Yang- Mills (SYM) theory [42|. In a 
constant medium, this approach leads to an depen- 
dence of mean energy loss. In the following, we will refer 
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to this approach by the label 'AdS'. 

We have refrained from testing a third class of models 
of energy loss by elastic scattering of a hard parton with 
medium constituents, which, as such scatterings arc in- 
coherent, would result in a linear dependence on the path 
length X in a constant medium. However, a large contri- 
bution of such processes to the total ener gy l oss can be 
ruled out already by the data for jislfiij. 

A key quantity in both models is the quenching weight, 
i.e. the energy loss probability distribution P{AE) given 
the path of a parton through the medium. In both models 
this is obtained by calculating the integrated virtuality 
transfer from the medium to the hard parton and 
the characteristic medium-induced gluon energy Wc by 
line integrals along the hard parton trajectory through 
the medium. Making use of a scaling law [45l|. we use 
Qs and Wc with the numerical results of [41| to obtain 
P{AE\ujc,Qs) = P{AE)pi,th- What is different between 
the two models is the argument of the line integrals. 

In the ASW model the medium is characterized by 
a transport coefficient q which measures the ability of 
the medium to transfer virtuality per unit path length. 
We assume that this can be written as a function of the 
medium thermodynamic parameters and parton position 
along the trajectory ^ via the relation 

q{0 e3/4(^)(coshp(C) - sinh p(0 cos a(0) (3) 

between the local transport coefficient g(^) (specifying 
the quenching power of the medium), the energy density 
e and the local flow ra pid i ty p with angle a between flow 
and parton trajectory [43 143 . We view the parameter K 
as a tool to account for the uncertainty in the selection 
of the strong coupling as and possible non-perturbative 
effects increa sing the quenching power of the medium (see 
discussion in l48|) and adjust it such that the pionic Raa 
measured in central Au-Au collisions is reproduced. 

With this expression for we evaluate for each path 
through the medium (given by the initial vertex position 
ro = (xcj/o) in the transverse plane and the angle 
of the outgoing parton with the reaction plane) the line 
integrals 



and 



(4) 



(5) 



In the AdS model, these expressions are changed into [42j 
Qs{roA) = K j d^^T\0 (6) 



and 



c^e(ro,0) = i^ j d^e T\0, 



(7) 



where T(^) the local temperature of the medium and 
K is again a (different) free parameter, to be adjusted 
such that the pionic Raa for central Au-Au collisions is 
reproduced. 

From the energy loss distribution for a given single 
path we can define the averaged energy loss probability 
distribution for a given angle (j) as 



oo roc 



— oo — oo 



{P{AE))^= dxo dyoP{xo,yo)P{AE)p,ti., (8) 



where for given impact parameter b the probability den- 
sity of hard vertices in the transverse plane P{xo,yo) is 
given by the product of the nuclear profile functions as 

p, ^ TA{r,+b/2)TA{ro-b/2) 
P[xo,yo) = tt; — tt: > (9) 

-lAA{b) 

and the nuclear thickness function is given in terms of 
the Woods-Saxon nuclear density PA{f, z) as 



TAir) 



dz pA{r,z). 



(10) 



We calculate the momentum spectrum of hard partons 
in leading order perturbative QCD (LO pQCD) (explicit 
expressions are given in j48| and references therein). The 
medium-modified perturbative production of hadrons at 
angle cj) can then be computed from the expression 



da. 



AA^h+X 
Tiled ^ 



E 



^ AA^f+X 

— ®(p(Ai?))^®i?];-,(z, pI), 

(11) 

with DJ°^f^(z, Hp) the fragmentation function with mo- 
mentum fraction z at scale fi^ (49j . From this we com- 
pute the nuclear modification function Raa vs. reaction 
plane as 



Raa{Pt,V,' 



dN\J dPrdydcj) 
TAA{b) daPP/dPrdydcj)' 



(12) 



The suppression of back-to-back high-p^ hadron cor- 
relations is computed in a Monte-Carlo (MC) framework 
[50l |. We start from the expression for the production of 
two hard partons k, I in LO pQCD which is described by 



da 



AB^kl+X 



dp\dyidy2 



= 2_,^i-^*M(^i'<3^)^2/j7i3(a::2,Q^)- 

ij 



dt 



(13) 

where A and B stand for the colliding objects (protons 
or nuclei) and yi{2) is the rapidity of parton fc(/). The 
distribution function of a parton type i in A at a mo- 
mentum fraction xi and a factorization scale Q pt is 
fi/A{xi,Q'^)- The distribution functions are different for 
free protons [5l|, [s^] and nucleons in nuclei [H, [13] . The 
fractional momenta of the colliding partons i, j are given 
by xi.2 = ^ (exp[±yi] + exp[±y2])- 

By sampling this expression, we generate events of 
back-to-back parton pairs which are placed on a vertex 
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sampled according to Eq. ([9|) for given orientation cf) with 
respect to the reaction plane. Given (rg, </>), we compute 
P{AE) for both partons according to the procedure out- 
lined above and sample the distribution to obtain the 
energy loss for the given event. 

Finally, we convert the simulated partons into hadrons. 
Note that this cannot be done using a fragmenta- 
tion function as in Eq. PT|) since takes a 
hadronic energy scale /xf as argument and measures the 
inclusive hadron yield, whereas we are interested in the 
yield of leading hadrons given a partonic energy scale. 

More precisely, in order to determine if there is a trig- 
ger hadron above a given threshold, given a parton k 
with momentum pT, we need to sample A^~^^{zi,pt), 
i.e. the probability distribution to find a hadron h from 
the parton k where h is the most energetic hadron of the 
shower and carries the momentum Pt = zi -pT- In the 
following, we make the assumption that the hadroniza- 
tion process itself, at least for the leading hadrons of a 
shower, happens well outside the medium. As a con- 
sequence, we neglect any interaction of formed hadrons 
with the medium. The time scale for hadronization of 
a hadron h in its rest frame can be estimated by the 
inverse hadron mass, Th ~ l/rrih', boosting this expres- 
sion to the lab frame one finds Tt ~ Eh/m\. Inserting 
a hard scale of 6 GeV or more for the hadron energy 
and the pion mass in the denominator (as pions con- 
stitute the bulk of hadron production), this assumption 
seems well justified. We extract Ai{zi,pt) and the con- 
ditional probability ^2(^1, Z2,pt) to find the second most 
energetic hadron at momentum fraction Z2 given that the 
most energetic hadron was found with fraction zi from 
HERWIG i53] . After hadronization, we check if the most 
energetic hadron fulfills a given trigger condition and, 
if yes, we count the yield in various momentum bins of 
hadrons back-to-back with the trigger. Finally, we obtain 
the suppression factor Iaa{4') for given trigger and asso- 
ciate momentum windows by dividing by the per-trigger 
yields found with nucleon parton distributions [Sll, [s^ in 
the absence of a medium. The procedure is described in 
detail in 



IV. THE BULK FLUID MEDIUM EVOLUTION 
'SEEN' THROUGH HARD PROBES 

Hard partons undergoing energy loss do not probe the 
same properties of the bulk medium as soft hadrons, or 
they probe them in a different way. For example, while 
the coefficient V2 in the soft sector measures pressure gra- 
dients translating an initial spatial anisotropy in non- 
central collisions into a momentum-space anisotropy, the 
same coefficient for high Pt hadrons measures directly 
the spatial anisotropy through the difi'erent energy loss 
induced by different densities seen by hard partons as 
a function of their angle with the reaction plane. This 
difference in the underlying physics is the reason why 
we prefer to present and discuss our results in terms of 



Raa (</') rather than in terms of the mean Raa and V2 at 
high Pt for a centrality class. While both choices contain 
the same information, we feel that Raa{<P) emphasizes 
the underlying suppression process. 

To give a second example, while mT^-spcctra of soft 
hadrons are rather sensitive to the late-time hadronic 
evolution of the medium and the amount of flow created 
during the hadronic evolution, the medium modification 
of hard probes is not at all sensitive to late time dynam- 
ics. The reason is that hard partons propagate through 
the medium with the speed of light, and thus typically 
escape from the medium at time scales of order of the size 
of the overlap region. This is especially true for observed 
hadrons, which have a bias to be produced relatively close 
to the surface [ist . 

In the following, we discuss some features of hydrody- 
namical models that are likely candidates to be probed 
by energy loss. 



A. Medium properties potentially probed by 
energy loss 

As apparent from Eqs. ([U to ([7]), the medium is probed 
by the energy loss models through path-length-weighted 
line integrals over the medium energy density e or tem- 
perature T . The value of these integrals thus depends 
on 

• the lower limit of the integral, corresponding to 
the time at which secondary particle production 
starts to be important enough to induce energy 
loss. Usually, the equilibration time tq of the hy- 
drodynamical model is used here, but there is no 
reason in principle why a non-equilibrated medium, 
if sufficiently dense, could not induce energy loss. 
However, the sensitivity to the lower limit is ex- 
pected to be comparatively weak, as factors of ^ 
or in the Wc integrals suppress this region. In 
a Bjorken model, q would diverge as 1/^ for small 
times, thus cancelling a factor ^ of the suppres- 
sion. However, prior to equilibration the density of 
particles in the medium off which the hard parton 
can scatter, thereby inducing it to radiate gluons, 
must generically be smaller than in thermal equilib- 
rium (which maximizes entropy and is thus a state 
of maximum particle density for given energy den- 
sity), so this cancellation cannot be perfect, and we 
expect (within reasonable limits) a weak sensitivity 
to the choice of the initial time. 

• the upper limit of the integral, corresponding to the 
time scale at which the hard parton is no longer sur- 
rounded by a medium. This is usually assumed to 
be given by the location of the Cooper-Frye sur- 
face in a hydrodynamical model beyond which the 
medium is no longer coupled but free-streaming. 
The Coopcr-Fryc surface in turn is often defined 
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to be an isothermal surface. Since the Cooper- 
Frye prescription is clearly an idealization of the 
real physics of the system boundary, there is again 
no strong reason to identify the upper integration 
boundary with this surface, as there could be en- 
ergy loss of a parton in a weakly coupled hadronic 
halo. However, again in practice the sensitivity to 
the detailed choice of the parameter is parametri- 
cally weak. While the factors ^ or tend to en- 
hance late-time contributions to the integral, the 
medium density at late times or large distances 
from the center dilutes eventually like due 
to both longitudinal and transverse flow. Thus, 
beyond a point we do not expect that our results 
depend strongly on the medium boundary choice. 

• the functional form of the integrand itself. Note 
that what is probed by the parton is not the 
medium density at any given proper time, but 
rather the medium density along the light cone. 
However, since interference effects suppress early 
time energy loss and the onset of transverse flow 
suppresses late time energy loss, effectively the 
mean energy loss per unit time dE/dx reaches a rel- 
atively sharp maximum around Tpcak ~ 3 — 4 fm/c 
(cf. j56|'). Thus, to first approximation, energy 
loss probes the density distribution of the medium 
around Tpeak- In particular, this implies that there 
is a finite time for processes like viscous entropy 
production or decay of the spatial deformation by 
pressure-generated anisotropic flow to modify the 
state of the system before it is probed by energy 
loss. 

• finally, when comparing the distribution ^ of pri- 
mary production vertices with any transverse dis- 
tribution of matter in a hydrodynamical model, one 
will find that there is always a non-zero probability 
to find a hard vertex outside the medium and hence 
partons which never experience energy loss. Since 
such 'halo-partons' never probe the medium, they 
show no correlation with the reaction plane angle 
and hence their presence will dilute any dependence 
of Raa on 0. This effect scales with the size of the 
hydrodynamical medium, i.e. with the assumed ini- 
tial extent of the Cooper-Frye surface. For very 
large media this 'halo effect' is suppressed, but for 
drastic assumptions (such as energy loss only in the 
QGP) it is not negligible. 

In the following, wc will discuss these points in more 
detail in the context of different hydrodynamical models. 

B. Viscous entropy production 

In both ideal and viscous fiuid dynamics, the entropy 
of the final state is fixed by the observed multiplicity and 
chemical composition of the emitted hadrons. But only 



in ideal fluid dynamics this final entropy agrees with its 
initial value. In viscous hydrodynamics, viscous heat- 
ing causes the entropy to increase with time, which im- 
plies that, for the same final multiplicity, the viscous fluid 
starts at lower entropy and parton density than the ideal 
fluid. The rate of viscous entropy production grows with 
the expansion rate of the flreball. A boost-invariant lon- 
gitudinal expansion profile, as initially assumed in all ver- 
sions of the hydrodyamic model studied here, leads to an 
expansion rate that diverges like l/r at early times r. 
Thus, most of the viscous entropy production happens 
near the beginning of the expansion (T^]. Viscous heat- 
ing causes the fireball density to decrease initially more 
slowly with time than for an ideal fluid. This has im- 
plications for the density profile seen by a hard parton 
propagating on the lightconc. 

Wc illustrate this effect in Fig. U where we show the 
integrand of Eq. ([S]), ^ q{£,), as a function of ^ for a parton 
propagating in-planc from the medium center in a hydro- 
dynamical background corresponding to Au-Au collisions 
at impact parameter b ~ 7.49 fm, for four different hy- 
drodynamical models. 




^[fm] 

FIG. 4: (Color online) The integrand of Eq.([5| as a function 
of parton distance ^ along the light-cone for four different 
hydrodynamical models, shown for a parton propagating in- 
plane from the fireball center. 



While the ideal hydrodynamical models show a mono- 
tonical decrease for ^^(Oi the viscous models show an 
initial rise due to viscous entropy production, followed 
by a decrease due to flow-driven density dilution. This 
means that in viscous hydrodynamical models energy loss 
is generically somewhat shifted to later times. 

Note also that due to the lower entropy in the initial 
state in viscous hydrodynamics, if one wants to get the 
same Raa, the factor K in Eqs. ([3]), ([6|), ([7]) must be 
larger (by about a factor 2 for the models presented here) 
than in an ideal fluid model that leads to a similar flnal 
state. 
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C. Line integral limits 

As discussed above, there is no compelling reason why 
the hydrodynamical thermalization time should be equal 
to the initial time of energy loss, nor why the Cooper-Frye 
decoupling surface for soft hadrons should coincide with 
the boundary beyond which the hard parton no longer 
interacts with the medium. 

In order to test the sensitivity of our results to choices 
of these parameters different from the hydrodynamical 
values, without explicitly modeling the (rather compli- 
cated) dynamics prior to thermalization or the perturba- 
tive scattering of hard quarks or gluons with hadrons in 
the halo, we adopt the following procedure: We probe 
the response of the system to variations of the hydro- 
dynamical starting time tq and of the final temperature 
Tp, with the understanding that these parameters are not 
meant as physically reasonable choices within the hydro- 
dynamical framework, but rather serve to generate upper 
limits for the true conditions — certainly prior to ther- 
malization the actual energy density will be less than in 
the hydrodynamical extrapolation to early times. In that 
sense, the parameters tq and Tp are to be understood in 
this subsection as the initial energy loss time and the 
equivalent temperature corresponding to the edge of the 
medium that contributes to energy loss. 
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FIG. 5: (Color online) The ratio of RAA{<t> = ^) I Raa (0 = f ) 
for various combinations of the initial energy loss time tq and 
the equivalent temperature Tp of the medium edge, for the 
case of the (2-|-l)-d ideal fluid model. 

In Fig. [5] we show the results of a variation of tq and 
Tp on the observable spread between the in-plane and 
out-of-plane emission, for the case of the (2-f l)-d ideal 
hydrodynamical model. While the variation looks opti- 
cally significant, it is small on an absolute scale (note 
the suppressed zero!), at most 15% from the mean. This 
is less than, for example, the variation between different 
hydrodynamical models. This in itself is certainly re- 
assuring, as it quantifies the uncertainty in choosing the 
proper line integral limits. It is also readily apparent that 
the spread is typically maximized for 'reasonable' choices 
of the freeze-out temperature, thus there is no evidence 



for a need to choose dramatically different last-scattering 
surfaces for soft and hard particles. 

Note that there is a systematic trend that large tq (i.e. 
delayed parton energy loss) leads to an increased in-plane 
vs. out-of-plane spread. This fits well into a pattern that 
shifting the strength of the mean energy loss per unit 
length dE/dx to later times leads to an increased ratio 
of out-of-plane vs. in-plane suppression. We will return 
to this issue in more detail later. 



V. RESULTS 

In order to illustrate how the differences between the 
media computed with different hydrodynamical models 
are probed by energy loss, we show results for Raa for 
in-plane and out-of-plane emission for various collisions 
centralities, using both the ASW and AdS energy loss 
frameworks and comparing with PHENIX data [531 ■ 
all calculations, a single free parameter K has been ad- 
justed such that a good description of Raa in 0-10% 
most central collisions is achieved. We therefore refrain 
from showing any results for central collisions, as they 
are virtually identical for all models [13] . 

In Fig.|n]we show results for the ASW energy loss mod- 
els from mid-peripheral collisions in the 20-30%, 30-40%, 
40-50% and 50-60% centrality classes. While all hydrody- 
namical models reproduce well the centrality dependence 
of the angular averaged (mean) Raa value, showing very 
weak Py-dependence, the spread between in-plane and 
out-of-plane emission appears generally too small (with 
the possible exception of the (3-|-l)-d ideal fiuid model). 
This is especially apparent in the 40-50% centrality class. 
Note that the successful description of the centrality de- 
pendence of the mean Raa is not trivial, since this probes 
the path length dependence of energy loss: in an elastic 
energy loss model with linear L dependence, even the 
mean Raa does not correctly extrapolate from central to 
peripheral collisions [4^ . 

Fig. [7] shows a similar comparison as Fig. [5] for the 
strong coupling AdS energy loss model which features 
a stronger L dependence than the ASW model (~ 
instead of ~ L^). Here there is a tendency for the mean 
Raa to lie above the more peripheral data, but the spread 
between in-plane and out-of-plane emission is generically 
larger and agrees better with the data. 

Let us comment on some trends: In all cases, the 
in-plane emission shows a stronger dependence on the 
underlying hydrodynamical model than the out-of-plane 
emission. This is to be expected, as the conditions for 
emission out-of-plane (e.g. in terms of average in-medium 
path length) are always more similar to the conditions in 
central collisions where in all models K has been adjusted 
to describe the data, whereas the conditions for emission 
in the reaction plane change more strongly with central- 
ity. 

The relatively small difference between fireballs result- 
ing from CGC and Glauber initial states is rather re- 
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FIG. 6: (Color online) The nuclear suppression factor Raa ss a function of Pt, shown in plane (solid) and out-of-plane (dashed). 
The calculation was done with the ASW perturbative radiative energy loss model for different hydrodynamical descriptions 
for the medium, describing 200 A GeV Au-Au collisions in four different centrality classes. PHENIX data [13 are shown for 
comparison. 



markablc. Given the stronger spatial anisotropy of the 
CGC initial state, one would expect to see this reflected 
in the anisotropy in high Pt parton energy loss. The so- 
lution to this puzzle seems to lie in the observation made 
earlier that energy loss probes the system at a timescale 
of ~ 3 — 4 fm, i.e. after the stronger pressure gradients 
of the CGC initial state already had some time to reduce 
the fireball eccentricity. 

With regard to the magnitude of the spread between 
in-plane and out-of-plane emission, we observe that con- 
sistently the (3-|-l)-d ideal hydro leads to the largest 
spread, followed by the viscous hydro models with CGC 
and Glauber initial conditions, while the (2-|-l)-d ideal 
fluid medium results in the smallest spread. We stress at 
this point that this observation cannot be directly linked 
to explicit modeling of the dynamics in z-dircction or to 
viscosity. We will discuss the causes for this ordering in 
the next section. 

In Fig. |8] we show the suppression factor Iaa for 
the away-side yield in triggered back-to-back correlations 
with a trigger momentum range of 4 — 7 GeV, using the 
ASW model. Note that the trigger range is not yet in 
a region where hadron production is dominated by the 
fragmentation of hard partons. It has been chosen to 



match with the range of an ongoing experimental analy- 
sis, keeping the mentioned caveat in mind. 

For 20-30% centrality, our statistics is not good enough 
to even cleanly separate in-plane from out-of-plane emis- 
sion. In the 50-60% centrality class, however, we essen- 
tially recover the ordering between models in the in-plane 
vs. out-of-plane spread observed before. 

Going to the AdS model in Fig. [9l the overall magni- 
tude of Iaa is somewhat different, but qualitatively the 
picture of the relative ordering of the spread between the 
different models remains unchanged. Given that com- 
putations of lAAi4') 8'i'6 rather involved and suffer from 
limited statistics, it is not clear that Iaa offers any real 
tomographic benefit over Raa- 

In order to illustrate more clearly how the energy loss 
probes the medium density evolution, we show in Figs. llOl 
and [TT] the conditional probability density to have the 
production vertex of the hard parton at position {x,y), 
given that a hard hadron was observed. We show calcu- 
lations for the 20-30% centrality class using the ASW en- 
ergy loss model; the qualitative features and differences 
between different hydrodynamical models seen in these 
plots repeat themselves in other centrality classes and are 
slightly enhanced for the AdS energy loss model which 
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FIG. 7: (Color online) Same as Fig. [6]except that the calculation was done with the AdS strong coupling radiative energy loss 
model. 



wc do not show here. 

Ill the absence of a medium, this production vertex 
density distribution is given by the binary collision dis- 
tribution Eq. ([9]). The medium biases the distribution 
in a characteristic way towards the surface, as partons 
produced in the dense medium core are unlilcely to es- 
cape with a significant fraction of their momentum left. 
The plots have been obtained by binning the distribution 
of triggered events in the MC code for the computation 
of back-to-back hadron correlations, but the same condi- 
tional P{x, y) for the trigger distribution is also underly- 
ing the single hadron suppression Raa- 

As expected, the degree of surface bias is much larger 
for out-of-plane than for in-plane emission, reflecting the 
larger degree of suppression seen out-of-plane. Compar- 
ing the y-position of the maximal out-of-plane emissiv- 
ity, a stronger degree of surface bias in the (2-|-l)-d ideal 
vs. the (3-t-l)-d ideal hydrodynamics is readily appar- 
ent. The two viscous models reflect their different den- 
sity evolutions in a less straightforward way, for example 
in the different shape of the emissivity maximum or in 
the distortion of the outer contour lines. As previously 
observed in [l^ , there is no evidence for strictly surface- 
biased emission or a corresponding tangential bias. In all 
models, a significant fraction of observed hadrons origi- 
nates from the medium core. 



VI. INTERPRETATION OF THE RESULTS 

Arguably the most important question in light of the 
magnitude of the measured spread between in-plane and 
out-of-plane emission is whether the data require a strong 
coupling description of energy loss along the lines of 
the AdS model or if, given a suitable description of the 
medium, perturbative QCD is able to account for the 
data. A related issue of similar importance is what con- 
straints for hydrodynamical models can be derived from 
measurements of hard probes that cannot also be gained 
from bulk matter data. 

Let us discuss these questions in view of our findings. 
It is a fortunate accident that the Cooper-Frye surfaces 
(see Fig. [3]) of all (2-|-l)-d models studied here are almost 
identical. Since we can easily test how RAAi4') is changed 
when we make tq (here denoting the starting time for en- 
ergy loss) equal in these runs (cf. Fig. [5]), we have prac- 
tically eliminated any effect of different in-medium path 
lengths and can ascribe any differences in the resulting 
suppression patterns to the density distributions probed 
along the hard parton path. 

We note that the (2-|-l)-d ideal and vGlb hydrodynam- 
ics both share (to good approximation) the same shape of 
the initial state whereas the vCGC hydrodynamics has 
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FIG. 8: (Color online) The suppression factor Iaa of the away-side per-trigger yield, calculated in the ASW perturbative 
radiative energy loss model for 4-7 GeV trigger momentum and shown as a function of the away-side momentum Pt- Squares 
show Iaa{(I> = 0) (in- plane), circles show Iaa (0= f) (out-of-plane) . Calculations are done for four hydrodynamical models for 
200 A GeV Au- Au collisions for two different centralities. 
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FIG. 9: (Color online) Same as Fig. [8] but for the AdS strong coupling radiative energy loss model. 



a different initial state. This enables us to disentangle 
different contributions to the in-plane vs. out-of-plane 
spread: about 50% of the difference in spread between the 
(2-|-l)-d ideal and vCGC hydrodynamics shown in Fig. [5] 
can be ascribed to the difference in the initial time tq 
(Fig. [S] indicates that a small tq is strongly disfavoured), 
about 35% of the difference results from viscosity, and 
the remaining 15% are caused by the difference between 
CGC and Glauber initial profiles. 

Turning to the (much larger) spread observed in the 
(3+l)-d hydrodynamics, we can surmise that most of the 
increase must arise from the different size of the freeze- 
out hypersurfacc (and the resulting different distribution 
of in- medium path lengths). This in turn implies that, 
once the parameter K has been adjusted to yield the 
same mean Raa as in the other models, the numerical 
value of q at each space-time point is much smaller in the 
(3-|-l)-d ideal hydro model than in the other models (see 
Fig.S]). Another way to state this is that in the (3-|-l)-d 
model partons travel on average a larger distance before 



they acquire the same amount of virtuality transfer from 
the medium. 

The common theme in all these findings is that the 
spread between in-plane and out-of-plane emission is in- 
creased whenever energy loss is shifted to later times. 
This may occur due to viscous heating during expansion, 
or in response to the choice of a large tq, or as the result 
of an path length weighting as in the AdS model — 
while the details differ, the net result is qualitatively the 
same in each of these cases. We note that this agrees 
qualitatively with the conclusions in [58| whose authors 
achieve a large in-plane vs. out-of-plane spread by requir- 
ing the hard parton to suffer the largest energy loss rate 
in a relatively thin shell of matter whose temperature is 
close to Tc- 

This observed connection between the time depen- 
dence of energy loss and the in-plane vs. out-of-plane 
spread requires an explanation, especially in light of 
the fact that elliptic flow decreases the initial spatial 
anisotropy of the system over time. The key to under- 
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FIG. 10: Probability density of finding a parton production vertex at {x, y) in the transverse plane, given an observed hard 
hadron with 4 GeV< Pt <7GeV. Calculations using the ASW energy loss model are shown for the (3+l)-d (top row) and 
(2+l)-d ideal hydrodynamical models (see text). In the case of in-plane emission (left panels) the hadron propagates to the 
—X direction and we use y -o- —y symmetrization for the plot. In the case of out-of-plane emission the hadron propagates to 
the -\-y direction and we used x o —x symmetrization. Countours are at linear intervals. 



standing this phenomenon lies in the realization that the 
spatial anisotropy exists on a scale of the order of the 
spatial size of the transverse overlap region, i.e. several 
fm. It is thus a global property of the medium. If energy 
loss were significant only at very early times t ^ Ifm/c 
(as approximately true for clastic energy loss), it could 
not resolve any phenomenon on a distance scale d S> 1 
fm and would thus probe the medium properties only 
locally. In other words, in such a model most high pT 
partons would be blind to the spatial anisotropy (see dis- 
cussion in [i^), and consequently the observed spread 



would be very small. 

While the above argument is strictly true only for a 
homogeneous medium, it still holds qualitatively for an 
inhomogeneous and hydrodynamically evolving medium; 
The spread between in-plane and out-of-plane emission 
increases not because energy loss is shifted to late times, 
but because the typical time scale over which energy 
loss is strong approaches the global spatial scales of the 
medium, and hence the partons are increasingly able to 
probe the medium globally. 

Unfortunately this finding does not translate into 



14 



2+ld viscous CGC, 20-30%, in plane 2+ld viscous CGC, 20-30%, out of plane 

6i 1 . . . — T . . . . 1 1 6i 





2+ld viscous Glauber, 20-30%, in plane 2+ld viscous Glauber, 20-30%, out of plane 





FIG. 11: Same as Fig. 1101 but for the viscous hydrodynamical models with CGC-fKLN (top row) and Glauber (bottom row) 
initial conditions (see text). 



straightforward constraints for the medium evolution 
model. Rather than being dominated by a single large 
effect, RAAi'P) appears to be sensitive to a number of ef- 
fects of roughly equal importance. However, while some 
constraints arc qualitatively similar to what was found in 
the soft sector (for example, a CGC initial profile causes 
both a larger V2 in the soft sector and a larger azimuthal 
variation of Raa in the hard sector), others are qualita- 
tively different and potentially more valuable (viscosity 
causes smaller V2 for the bulk matter but a larger az- 
imuthal variation of Raa^^I^))- 

Given that our results do not exhaust the parameter 
space of hydrodynamical evolutions compatible with bulk 
data, it is entirely conceivable that an evolution can be 



found for which both magnitude and spread of RAA{(t') 
are described within pQCD dynamics. Presumably, a 
viscous (3-f l)-d hydrodynamics with a CGC initial state, 
late thermalization and low freeze-out temperature would 
be a good candidate. At this point, we sec no reason to 
conclude that the data require a strong-coupling model. 

VII. SUMMARY 

The goal of this paper was to develop a better un- 
derstanding of the tomographic power of parton energy 
loss measurements in relativistic heavy-ion collisions as 
probes of the medium created in these collisions and/or 
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of the mechanism by which hard partons lose energy in 
such a medium. To this end we presented a study of 
the dependence on the angle with the reaction plane (p 
of the nuclear suppression ratio Raa{<P) and the away- 
side per-trigger yield /yiA(</>) in triggered back-to-back 
correlations, for four different hydrodynamic models of 
the fireball evolution and two different models for the 
path length dependence of parton energy loss. The mod- 
els were tightly constrained by ensuring that parameters 
were chosen such that the soft-hadron transverse momen- 
tum spectra are well described for all collision centrali- 
ties, and that all models yield the same nuclear suppres- 
sion factor Raa for pious in central Au-Au collisions. 
We then studied differences in the dependence of Raa 
and Iaa on collision centrality, transverse momentum 
Pt, and azimuthal emission angle 4> relative to the re- 
action plane. 

We found that, after tuning the models to reproduce 
the correct pion Raa in central collisions, they all gave 
approximately identical results for the collision central- 
ity dependence of the azimuthally averaged Raa{Pt), 
featuring weak Pr-dependence and yielding good qual- 
itative agreement with experimental data. The tomo- 
graphic power of this azimuthally averaged quantity for 
distinguishing between different (realistic) medium mod- 
els is therefore low. A stronger path length dependence 
of the energy loss instead of ~ i^) produces a 

slightly stronger impact parameter dependence of the 0- 
averaged Raa , but small differences between the different 
hydodynamic models for the medium interfere with this 
tendency, making it difficult to disentangle these effects. 
The centrality dependence of the azimuthally averaged 
Iaa is a priori a bit more promising since it was found to 
react more strongly to changes in the path length depen- 
dence of the energy loss. On the other hand, to measure 
this quantity with good statistical precision is much more 
difficult. Furthermore, it is not obvious from our studies 
that such a measurement will contribute any useful in- 
formation about the medium that would help distinguish 
between different hydrodynamic evolution models. The 
azimuthally averaged Raa is found to provide practically 
no such discriminating power. 

The measurement of the in-plane vs. out-of-plane vari- 
ation of Raa{4') provides much better discrimination 
[65j . We confirm earlier findings that it is not easy to 
reproduce the relatively large oscillation amplitude of 
i?AA(</') found by the PHENIX experiment. We iden- 
tified several mechanisms that help to increase the in- 
plane vs. out-of-plane spread of Raa and, in combina- 
tion, may be able to explain the data: (i) A stronger 
path length dependence of parton energy loss, combined 
with (ii) a delayed beginning of the energy loss action, 
caused by the need for allowing the scattering centers 
in the medium that induce the hard parton to radiate 



energy to decohere from the initial state wave function 
of the colliding nuclei, (iii) a delay of the flow-induced 
dilution of the medium density by viscous heating, and 
(iv) a larger fireball eccentricity by intializing the hy- 
drodynamic evolution with CGC-fKLN initial conditions 
rather than the less deformed Glauber model initial den- 
sity profile. We found that the combination of effects (i) 
and (ii) accounts for about 50% of the azimuthal spread 
found in our calculations, viscous heating contributes an- 
other 35% of the effect, with the remaining 15% arising 
from different initial fireball eccentricities. 

Analyzing the reasons why the mechanisms (i)-(iii) 
cause a larger in-plane vs. out-of-plane spread of Raa, 
we found that they all shift the weight of the energy 
loss towards later times along the path of the parton. 
This allows the parton energy loss to probe the spatial 
anisotropy of the medium (which is a global fireball prop- 
erty) on a global length scale, rather than probing the 
properties of the medium locally in the vicinity of the 
production vertex only. If the global spatial anisotropy 
of the medium can be probed by a parton, the connection 
between spatial medium anisotropy and final state hard 
parton momentum anisotropy is strongest. 

We would like to point specifically to the role of vis- 
cosity in this context: We are beginning to see signif- 
icant roles played by the (necessarily non-zero) shear 
viscosity of the quark-gluon plasma even though its ab- 
solute value (expressed through the dimensionless ratio 
^ = 0(1— 3) X ^) is almost as small as theoretically 
possible In [2lj we found that ^ « (2—3) x ^ gives 
the best fit to the pion and proton momentum spectra 
in 200AGeV Au-Au collisions, and here we found that 
a similar value contributes significantly to the in-plane 
vs. out-of-plane spread of Raa for pious. We hope and 
expect that a combined analysis of all RHIC data on soft 
hadron production and hard parton medium modifica- 
tion will eventually lead to an accurate determination of 
the quark-gluon plasma viscosity. 
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